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ABSTRACT 

We present photometry of six transits of the exoplanet X0-2b. By combin- 
ing the light-curve analysis with theoretical isochrones to determine the stellar 
properties, we find the planetary radius to be 0.996^^^ i?j up and the planetary 
mass to be 0.565 ± 0.054 M Jup . These results are consistent with those reported 
previously, and are also consistent with theoretical models for gas giant planets. 
The mid-transit times are accurate to within 1 min and are consistent with a con- 
stant period. However, the period we derive differs by 2.5a from the previously 
published period. More data are needed to tell whether the period is actually 
variable (as it would be in the presence of an additional body) or if the timing 
errors have been underestimated. 

Subject headings: planetary systems — stars: individual (XO-2, GSC 03413-00005) 



1. Introduction 

Observations of exoplanetary transits have provided the first empirical information 
about the internal structure, composit ion, surface temperature and atmospheric dynam- 



ics of planets outside the Solar system (ICharbonneau et al.l 120071 ). It is hoped that precise 
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measurements of transit t imes and durations will also provide a new channel for the detec- 
tion of low- mass planets (IHolman fc Murray! 120051 ; lAgol et al.ll2005l ). Among the immense 
amount of information these systems may deliver, the most fundamental parameter is the 
radius of the planet. Accurate radius determination is not a trivial task, and in some cases 
after a transiting planet has been announced, more precise data have led to significant revi- 
sions on their radius, with consequences for the ories of planetary interiors and atmospheres 
(see, e.g. Iwinn et aliboosl : I Johnson et aDl2008h . 



The Transit Light Curve (TLC) project is an effort to gather precise photometry of 
exoplanetary transits, in order to determine the fundamental properties of these planets and 
their parent stars as accurately as possible, and to seek evidence of undetected planets in 
the pattern of transit times and durations. 



The present paper is concerned with the X O-2 system (IBurke et al 



the second system discovered by the XO project (IMcCullough et al.ll2006l ). The XO-2 star 



2007|). This was 



is a nearby Solar-type star with a binary companion (identified through its common proper 
motion), and has a high metallicity ([Fe/H] = 0.45) compared to other planet-hosting stars. 
The XO-2b planet is a normal hot Jupiter (if the word "normal" may be applied to such 
interesting obj ects), with an orbi tal period P = 2.615857 d, mass M = 0.57M Jup , and radius 
R = 0.98R 3np faurke et al. l l2007h . 



In this work we present and analyze differential photometry of six transits of XO-2b. 
We measured the mid-transit time for each of the six events, and by combining the data, we 
obtained new and independent estimates for the system parameters. This paper is organized 
as follows. In § [2] we describe the observations and data reduction procedure. In § [3] we 
describe the model and the techniques we used to estimate the physical parameters of this 
system. In § 0] we present the results for the planetary parameters and mid-transit times. 
Finally, in § we present a brief summary of this work. 



2. Observations and Data Reduction 

We observed four complete transits, and one partial transit, with the 1.2 m (48-inch) 
telescope at the Fred L. Whipple Observatory (FLWO) on Mount Hopkins, in Arizona. The 
FLWO 1.2m telescope is equipped with KeplerCam, a monolithic 4Kx4K Fairchild 486 CCD 
that gives a 23' x 23' field and a pixel size of 0.68" in 2x2 binned mode. We used a Sloan z 
band filter to minimize the effects of limb darkening on the shape of the transit light curve. 

We observed one complete transit with the Centurion 0.46 m (18-inch) telescope, at 
the Wise Observatory, in Israel. The Centurion telescope is equipped with a Santa Barbara 
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Instrument Group (SBIG) ST-10 XME USB CCD camera. This thermoelectrically-cooled 
chip has 2184 x 1472 pixels each 6.8 jum wide, which convert to 1.1 arcsec at the //2.8 focus 
of the telescope. The chip offers a 40.5' x 27.3' field of view. Observations from Wise were 
made with no filter, because this camera has no filter wheel. The 18-inc h telescope CCD h as 
an overall response similar to a "wide-R" band (for further details, see lBrosch et al.ll2008l ). 

Relative aperture photometry was performed for XO-2 and 13 nearby comparison stars, 
including the binary companion of XO-2 (7 stars were used for the Wise data set). The choice 
of comparison stars was iterative; we began with a longer list, and those stars that showed 
unusual noise or variability were removed. The 13 comparison stars range in brightness from 
50-200% of the brightness of XO-2. The sum of the fluxes of the comparison stars was taken 
to be the reference signal, and the flux of XO-2 was divided by this reference signal. The 
root-mean-square (rms) noise in the relative flux of XO-2 among the six light curves ranges 
between 0.0013 and 0.0036. Table [1] gives a summary of the six photometric observations, 
with information about the dates, epochs, band, exposure times and RMS. Figure [1] shows 
the observed light curves. A differential airmass correction was applied to the plotted light 
curves. This correction was determined as part of the fitting process that is described in the 
next section. Table [2] has an abbreviated versions of one of the light curves. We intend for 
the complete photometric data set to appear in the electronic version of the journal. 



3. Data Modeling 
3.1. Light Curve Analysis 



We modeled the light curves using the analytic formulas of iMandel fc Agoll (120021 ). The 
orbital period was held fixed using the value of Burke et al. (2007) and the planetary orbit was 
assumed to be circular. The fitted parameters were the radius ratio R p /R+, the normalized 
semi-major axis a/R+, the impact parameter b, the mid-transit time T c , and two parameters 
that take into account the zero-point and slope of the differential airmass correction. We 
assumed a quadratic limb-darkening law, with coefficients u\ and u<i- Through a trial fit, we 
found that the data are not precise enough to fit for both limb-darkening coefficients. The 
linear combination 2u\ + «2 is well-constrained, and the orthogonal combination u\ — 2m 2 is 
poorly constrained. For this reason we allowed 2u\ + u-2 to be a free par ameter, and held 
ui — 2w 2 fixed. We fixed u\ — 2w 2 at the value tabu lated by Claret ( 200 4)^1. for a star with 
the metallicity, gravity, and temperature derived by lBurke et al.l (120071 ). 



^romlClard (120041 ). u x = 0.292 and u 2 = 0.299 
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The fitting statistic was 



Xlc 



E 



yobs jmod 



(J; 



where f° bs and /j 71 are the observed and modeled relative fluxes observed at time j, and aj 
is the uncertainty of the data points as described below. To set the appropriate data weights, 
we used a procedure that attempts to account for any time-correlated ("red") noise, at least 
approximately. For each light curve we found an initial best-fitting model (using the out-of- 
transit RMS as the flux weight) and calculated <Ti, the standard deviation of the unbinned 
residuals between the observed and calculated fluxes. Next we averaged the residuals into 
M bins of N points and calculated the standard deviation of the binned residuals. In the 
absence of red noise, we would have expected 



0"! 



M 



Vnv m-i 



(2) 



but often the measured value of is larger than this by a factor (3 due to time-correlated 
noise. In such cases the number of effectively independent data points is smaller than the 
actual number of data points. In Eq. [I] we used Oj — <j\ x /3. We found that (3 depends 
weakly on the choice of averaging time r, generally rising to an asymptotic value at r ~ 10 
minutes. The value (3 was close to unity except for the light curves of UT 2008 Feb 25 
(0 = 1.2) and UT 2008 Mar 06 (fi = 1.5). 

The best-fit parameters and their uncertainti es were obt aine d using a Markov Ch ain 
Monte Carlo (MCMC) procedure. As described by iFordl (120051 ) and lHolman et al.l (120061 ). in 
this method a random process is used to create a sequence of points in the parameter space 
that approximates the studied probability distribution. This sequence or chain is generated 
by a jump function that adds a Gaussian random number to each parameter. The jump 
is executed if the new point has a xf c l° wer than the previous point. If xf c * s l ar g er ; the 
jump is made with a probability equal to exp(— Axf c /2). If the jump is not made, the new 
point is a copy of the previous one. The relative sizes of the perturbation were set equal to 
the approximate la uncertain t ies ob tained by direct inspection of xf c across the parameter 
space, as done in lBeatty et al.l (120071 ). The sizes of the jumps are set by requiring that ~25% 
of the jumps are accepted. As a first step, four independent chains of 2 x 10 5 points were 
created for each light-curve, discarding the first 20 % of points to minimize transient effects. 
The four chains were then superposed to create one sequence of points. The Gelman & 
Rubin (1992) R statistic was always within 0.2% of unity for each parameter, a sign of good 
mixing and convergence. The cumulative posterior probability is a normalized histogram of 
the MCMC sample values. To derive the confidence interval for a parameter, the MCMC 
sets were sorted by the parameter of interest and then we determined the 15.85%, 50%, and 
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84.15% points of the cumulative posterior distribution. The 50% point (i.e. the median) was 
taken to be the "best-fitting value" and the interval between the 15.85% and 84.15% points 
was taken to be the 68.3% confidence interval. 

The final values for the transit parameters R p /R+, a/R+ and b were obtained by analyz- 
ing the combined data from the five light curves observed with FLWO /Kepler Cam, following 
the same procedure used for the individual transits. A 2-minute binned version of these com- 
bined light curves is shown in Figure [2J 

For the transit timings, we also performed two independent additional tests to check 
the accuracy of our results. For one of the tests we used a standard bootstrap simulation: 
we created 10 4 "realizations" of the data set by perturbing the best-fitting model with 
Gaussian random noise (with a standard deviation equal to the rms of the actual light 
curve), and minimizing xf c as a function of t c for each realization (Rp/R*, a/R* and b are not 
correlated with t c , and were therefore held fixed for these tests). The resulting collection of 
10 4 t imings for each tr ansit was taken to be the probability distribution of the timing (see, 
e.g., iPress et al.l 119921 ). For the second test we implemented the residual-permutation or 
"rosary-bead" method. For each light curve, we found the best-fitting model and computed 
the time series of N residuals. We then added these residuals to the model light curve 
after performing a cyclic permutation of the time indices. This is a variant of the bootstrap 
technique that preserves the temporal 
previously by many investigators (e.g., 



Bouchv et al. 


2005: 


„ 

Southworth 


2008) 



light curves under study, the timings obtained on the additional tests were consistent with 
our first set of results, and with similar error bars. The results of the light curve analysis 
are summarized in Table [31 



3.2. Determination of Absolute Dimensions 



As has been shown before by several authors, the only intrinsic properties of the star 
and planet that can be determine d directly from observed quantit ies on transiting systems 
are the mean density of the star (ISeager fc Mallen-Orne las 20 03|) and the surf ace gravity 



of the planet (ISozzetti et al.l 12007 



Holman et al. 



2007 



Southworth et al.l 12004 ) . In order 



to determine the individual masses and radii of the bodies, external information must be 
introduced. 

In th is case, we u s ed st ellar evolution models following the procedure and consider- 
(120081 ). For this purpose we rely on the spe ctroscopic t e mpera - 
5340 K) and metallicity ([Fe/H] = 0.45) obtained by burke et all J2007h . 



at ions of iTorres et al 
ture (T eS 
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but we adopt more conservative errors, no smaller than 0.05 dex in [Fe/H] and 80 K 
in T e g. The reason for this approach is to take into account the documented difficulty 
in obtaining accurat e values for stellar temperatures (see, eg.. Ifjlackwell &: Shallis 1977 ; 

2005h and metallicities 



Blackwel 



et al.lll980l; ICasagrande et al 



2006; 



Ramirez &; Melendez 



(see, eg. JFischer fc Valenti 20051 ; Gonzalez fc Laws! 20071 ; Santos et al.ll2004j ). Instead of us- 
ing the spectroscopic value for logg, we used the parameter ct/i?* that is closely related to 
the mean stellar density and it is provided by the transit light curve fit. The quantity a/R* 
can be obtained from the models using the following expression: 



a/R* 



G \ 1/3 P 2 / 3 



4vr 2 



(M, + M p ) 1/3 



(3) 



The mass of t he planet is not kno wn a priori, but a/R+ can be estimated using the value of 
M p derived by lBurke et al.l (120071 ). At the end of the procedure to be explained next, a new 
value of M p is obtained, which is used to repeat the process until convergence. 



The stellar evolut ion models used were those from the Yonsei-Yale series ( lYi et al.ll2001 



Demarque et al.l 12004 ) . These isochrones were interpolated to a fine grid in metallicity and 
age and compared point by point with the measured values of T e g and a/R+. Each point 
on the isochrones that was consistent with [Fe/H], T c g and a/R* within their errors, was 
recorded and a likelihood given by exp(— % 2 /2) was calculated, where 



A [Fe/H] 

C [Fe/H] 



+ 



AT, 



cff 



a T cff 



+ 



A(a/R* 



a/Rt, 



(4) 



In this expression the A symbols represent the difference between the observed and model 
values for each quantity. The best fit value for each stellar parameter was obtained by adding 
all matches, weighted by their corresponding normalized likelihood. We did not account for 
the varying density of stars on each isochrone, as the effect is generally small in the case 
of solar-type stars. The adopted errors for the fitted parameters (mass and age) come from 
their range among the accepted points on the isochrones. Figure [3] illustrates the location of 
the star in a diagram of a/R+ vs. effective temperature, similar to an H-R diagram, compared 
with isochrones from the Yale-Yonsei series. 

With M* known, we obtained M p and a by iteration of Newton's modified version of 
Kepler's third law and the mass function of the system: 



G 
4^2 



(M* + M p ) P 2 



M n 



(—) 

\2nG 



1/3 



sin i 



(M* + M p ) 



2/3 



(5) 
(6) 
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The value of a in combination with a/R* and R p /R* allowed us to obtain consistent values 
for R+ and R p . 

To estimate the errors for M p , .R* and R p , we used the MCMC chains generated in the 
course of modeling the transiting light curves. For each element of the chain a solution was 
calculated using random Gaussian values for P and K+ (using the observed values and their 
uncertainties as the center and standard deviation of the Gaussian distributions). In this 
way we obtained a probability distribution for M p , R+, and R p , from which we extracted 
the median and 68.3% confidence intervals and adopted them as best values and errors, 
respectively. 

It is important to note that our procedure places complete trust in the Yonsei-Yale 
isochrones. It assumes they are exactly correct. Any systematic errors in the isochrones are 
not accounted for in our error bars. As this star is very similar to the Sun, we expect these 
errors to be small. 



4. Results 

Table [3] gives all the measured planetary and stellar properties, together with the results 
for the light-curve parameters. The labels, explained in the caption, clarify which quantities 
were obtained from the literature, which were determined independently by our analysis, 
which quantities are functions of those independent parameters, and which quantities depend 
on our stellar isochrone analysis. 



4.1. Planetary and stellar Properties 



The light curve best-fit pa rameters a/R * : Rp/ R-k and b presented in this work are con- 
sistent with those obtained by iBurke et al.l (120071 ) in the discovery paper of XO-2b, with 
similar error bars for all of them. The precision of these parameters is dominated by the first 
and third light curves, which have the highest signal-to-noise ratio. For our stellar isochrone 
analysis, we used the values for T eS and [Fe/H] from the discovery paper about XO-2b, 
together with our value for a/R+. W e obtained a stellar mass value (M* = 0.971 ± 0.043) 
almost identical to the one found by IBurke et al.l (120071 ). In our case, the result does not 
depend on the assumed distance to the star, but only on T e s and [Fe/H]. We derived the 
planetary mass and radius plus the stellar radius directly from the stellar mass. None of these 
parameters show any significant discrepancy from those listed in the discovery paper. Our 
error bars are similar as well, but in this work we took into account the red noise in the light 
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curves and adopted larger errors for T e g and [Fe/H], and in this sense our error estimates are 
more conservative. In the case of the stellar radius .R*, the value obtained directly from the 
isochrone modeling and the one one derived from the stellar mass are identical, but the error 
bars on the latter are 25% smaller. Our st ellar surface gravity estimate (log g+ = 4.45 ±0.01) 
is independent from the one obtained by lBurke et al.l (120071 ). Both estimates of the surface 
gravity are in agreement, and our result based on a/R* is more precise. This illustrates the 
power of tran sit light curves to pin down the stellar mean density, as also exhibited for the 
XO-3 system Jwinn et aliboosh . 



As mentione d earlier, we determ ined the stellar properties following the procedure and 
considerations of iTorres et al.l (120081 ). therefore, we expected our results for XO-2b to be 
in good agreement with theirs. This is actually the case, given that our analysis of the 
transit light curves didn't deliver parameters significantly different to the ones obtained by 
Burke et all (120071 ). 



We f ound that the mass a nd radius for XO-2b are in good agreement with those pre- 
dicted by iFortney et al.l (120071 ) for a planet orbiting a solar-like star. Interpolation of the 
tabulated results presented on their work for the appropriate values of the XO-2 system 
(M p = 0.57Mj up , a = 0.037 AU) gives a theoretical planetary radius of 1.0i?j up . 



4.2. Transit Timings 

Using our 6 transit times, we computed an ephemeris independent of that of Burke et 
al. (2007), by fitting a linear function of transit epoch E, 

T C (E) = T + EP. (7) 

The result is T = 2, 454, 466.88514 ± 0.00019 [HJD] and P = 2.615819 ± 0.000014 days. 
The fit has x 2 = 5.76 with 4 degrees of freedom (iVd f = 4), indicating an acceptable fit (i.e. 
X 2 is within y/2N do{ of iV do f). However, the fitted period differs by 2.5a (3.3 s) from the 
value 2.615857 ± 0.000005 d that was determined by Burke et al. (2007). 

It is not clear how to interpret this discrepancy. It may be a genuine period variation 
produced by an additional body in the XO-2 system. However it seems at least as likely that 
the uncertainties in one or a few of the transit times have been underestimated. Our methods 
and cross-checks on the errors have already been described. As for Burke et al. (2007), 
one of their mid-transit times was based on observations with a 1.8m telescope, and has a 
quoted precision that is better than any of our 6 measurements (2, 454, 147.74902±0. 00020). 
The other 11 light curves came from observations with smaller telescopes. For those 11, 
Burke et al. (2007) estimated an error of 3 min, based on the scatter among 3 independent 
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measurements of one particular transit in March 2007. The relative quality of the individual 
light curves, which were gathered by different observers on different nights, was not taken 
into account. 

One might wonder if the discrepancy can be attributed to any particular data point. If 
we use only the most precise time of Burke et al. (2007) in combination with our 6 times, 
fitting a linear ephemeris gives \ 2 = 15.9 and A^of = 5, which is still unacceptable. If we 
omit the most precise time and instead use the other 11 times from Burke et al. (2007), 
adopting 3 min errors as they did, we find \ 2 = 20.9 and iVdof = 15, right on the margin 
of acceptability. If we use all of the times from Burke et al. (2007) and omit our last 2 
times, which were derived from the noisiest light curves, we find x 2 = 15.6 and iVdof = 14, 
a satisfactory fit. (Omitting either of our last two times individually does not result in a 
satisfactory fit.) No further conclusions can be drawn until more times have been measured, 
over a longer time baseline. 

For planning purposes, we recomputed the ephemeris based on all the available transit 
times (6 from this work, and 12 from Burke et al. 2007). The result is T = 2, 454, 466.88467± 
0.00013 [HJD] and P = 2.6158640 ± 0.0000016 d, with X 2 = 27.6 and N do{ = 16. Observers 
who are using these values to plan future observations may wish to inflate the error bars 
by a/ X 2 /Ndo{ = 1-31 to account for possible systematic effects due to either genuine period 
variations or underestimated measurement errors. 

A more detailed study of possible short-term transiting timing variations caused by 
additional orbiting bodies would requir e the o bservation of a larger number of conse cutive 



transits, as noted by lHolman fc Murray! (120051 ) and tested by iMiller-Ricci et al.l (120081 ). Any 



constraint derived from the present data would be of limited interest and beyond of the scope 
of this paper. 



5. Summary 

We have presented new photometric observations of six transits of the exoplanet XO-2b. 
The analysis of the data independently confirms the previously calculated properties of this 
planet, leaving little doubt about the orbit and geometric properties of the XO-2 system. 
We found that the mass and radius for XO-2b are in good agreement with standard models 
of gas giant planets. The six new transit timings yield a refined transit ephemeris and set 
an updated and precise reference for future searches of secondary bodies orbiting this star. 
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Table 1. Photometric Observations Summary 



UT Date 


Jan 01 2008 


Jan 14 2008 


Feb 12 2008 


Feb 25 2008 


Mar 04 2008 


Mar 06 2008 


Epoch 1 





5 


16 


21 


24 


25 


Observatory 


FLWO 


FLWO 


FLWO 


FLWO 


FLWO 


Wise 


Band 


Sloan 2 


Sloan 2 


Sloan 2 


Sloan 2 


Sloan 2 


clear 


Exposure [sec] 


30 


25 


25 


30 


25 


9 


FWHM median [arcscc] 


1.6 


1.5 


1.2 


1.4 


2.9 


1.3 


RMS [rcl. flux] 


0.0013 


0.0013 


0.0017 


0.0027 


0.0026 


0.0036 



Note. — (l)Numbcr of cycles elapsed since the initial transit studied in this work 
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Table 2. Transit Light Curve 



Jan 01 2008 UT 


HJD 


Relative Flux 


2454466.78551 


0.99807 


2454466.78602 


1.00079 


2454466.78654 


1.00180 


2454466.78705 


0.99955 



Note. — This is a sample entry 
of a full light curve. The complete 
versions are given on-line. 
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Table 3. Transit Timings and System Parameters for XO-2 



Parameter 




Value 


68.3% Conf. Limits 


Comment 


T c (0) 

c v / 


[HJD1 

L J 


2,454,466.88512 


0.00021 


A 


T c (5) 


[HJD] 


2,454,479.96393 


0.00039 


A 


T c (16) 


HJD 


2,454,508.73864 


0.00026 


A 


T c (21) 


[HJD] 


2,454,521.81778 


0.00072 


A 


T c (24) 


[HJD] 


2,454,529.66433 


0.00043 


A 


T c (25) 


HJD 

L 1 


2,454,532.27978 


0.00074 


A 


P 


[days] 


2.6158640 


0.0000016 


B 


To 


[HJD] 

L 1 


2,454,466.88467 


0.00012 


B 


-Rp/ -R* 




0.10485 


+0.00070, -0.00062 


A 


a/R+ 




8.13 


+0.09,-0.20 


A 


b 




0.16 


0.11 


A 


Ui 




0.250 


0.026 


A 


U2 




0.296 


0.013 


A 


Pstar 


[g cm- 3 ] 


1.484 


+0.051,-0.104 


C 


log g p 


[cm s~ 2 ] 


3.147 


+0.043, -0.048 


C 


M* 


[M Q ] 


0.971 


0.034 


D 


age,. 


[Gyr] 


6.3 


2.4 


D 


R* 


[Re] 


0.976 


+0.024, -0.016 


E 


log g* 


[cm seg~ 2 ] 


4.448 


+0.011,-0.021 


E 


M p 


[M Jup ] 


0.565 


0.054 


E 


R p 


[-Rjup] 


0.996 


+0.031,-0.018 


E 


a 


[AU] 


0.0368 


0.0004 


E 
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Note. - (A) obtained independen tly. (B) derived using transit timings 

from this w o rk an d from iBurke et al.l (120071 ). (C) calculated using K* from 



Burke et al.l (120071) .fD) deriv e d from isochrone modeling using values for T e g- 



and [Fe/H] from Burke et all (120071 ). (E) function of A, C and D 
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FLWO - 20080101UT 




o 



-3-2-10123 

time - T c (hrs) 

Fig. 1. — Individual transit light curves, after applying a correction for differential extinction. 
The FLWO transits were observed in the Sloan z band. No filter was used at Wise; the 
effective bandpass resembles wide R band. 
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Fig. 2. — A. Combined FLWO light curves. B. Best-fit model and 2 min binned combined 
light curve. The light curve observed from Wise is not included in this combined data set. 
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8000 7500 7000 6500 6000 5500 5000 4500 



T ef f I K ] 

Fig. 3. — Yonsei-Yale stellar isochrones for the modeling of the primary star, corresponding 
to ages of 0.6 Gyr, 1.6 Gyr, 2.6 Gyr, etc. (left to right), and the chemical composition 
indicated. The observed transit parameter a/R±, which is closely related to the mean stellar 
density, was used as an indicator of luminosity. The black dot and the filled box around it 
indicate the best values and uncertainties for a/R* and T eff . 
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Fig. 4. — Differences between observed and calculated transit times. On panel (a), the 6 TLC 
transit timings are shown. The time residuals are calculated from the ephemeris obtained 
using TLC transits only. Panel (b) shows the time residuals of all the transits reported on 
the discovery paper plus those presented in this work. The time residuals are calculated 
from the ephemeris obtained using all the available timings. 



